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We examine the different possible symmetries of the superconducting gap obtained by 
solving the Eliashberg equations. We consider an electron-phonon interaction in a strong 
coupling scenario. The Coulomb pseudopotential plays the crucial role of providing the 
repulsion needed to favour the d-wave symmetry. But the key parameter that allows very 
anisotropic solutions even with very strong coupling is the small angular range of the 
interaction due to predominantly electron-phonon forward scattering that is found in the 
high- Tc superconductors. We find both s and d-wave solutions whose stability depends 
mainly on the angular range of the interaction. 
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I. INTRODUCTION 

In the last few years, much experimental and theoretical effort has been spent in trying to understand 
the mechanism responsible for superconductivity in the high-Tc superconductors. Such a mechanism 
might be the BCS electron-phonon interaction, the exchange of paramagnons, or some Coulomb inter- 
action processes as in the t-J model or in the Pines pictureEl. The symmetry of the gap has been given 
much attention since it may clarify some aspects of the superconducting mechanism. It seems now quite 
established, at least in YBCO, that the symmetry of the gap is d-waveB. 

Our aim in this paper is to show that a standard electron-phonon attractive interaction in conjunction 
with a weak Coulomb repulsion can lead to a d-wave symmetry of the superconducting gap A. This 
goes against the common idea that an electron-phonon mechanism can only lead to an s-wave symmetry. 
Indeed, if the sign of the^interaction does not change, there^is a theorem in linear algebra saying that the 
gap will not have nodestl. Furthermore, it has been showno that a large electron-phonon coupling would 
wash out any anisotropy in the gap. It is true that in normal metals, the electron-phonon interaction is 
local due to the strong screening, thereby allowing averaging over the Fermi surface (FS)El. In the cuprates, 
the situation is different. The dielectric constant of the ionic background is very large for small frequencies, 
i.e. eion(w ^ 0) « 40 — 60 in YBCOQ (what we denote by eion is the dielectric constant obtained after 
substracting the Drude contribution of the conduction electrons at frequencies lower than the frequencies 
of the transverse optical modes; in this frequency range, e is nearly independent of to; this quantity is 
sometimes denoted CooQ). Since the response of the ions to an external charge is so strong, the response 
of the electrons, given normally by eei{q) = 1 + {nTp/qY with ni^p = Ane'^N{eF), will be decreased. 
Therefore, the dressed electronic dielectric constant ^'^J'^'^'^^^ ^^^^ -^^^ written as ^dressed f^^--^ = 1 -f- (gg; — 

l)/eion = ^ + K'^p/ici'^^ion)- The range of the interaction in /c-spacc for a uormal metal is givcn by kt_f (~ 
2kp), while in our case, ktf / \/^ion ~ 0.3kp. Therefore, the electron-phonon scattering is predominantly 
forward, i.e. the angular range of the interaction (for electrons at the FS) SO « ktf / {^ionkp) ~ tt/IO is 
very small. This electron-phonon interaction then allows for a very anisotropic s-wave gap and even a 
d-wave one, due to the repulsive Coulomb interaction (since the sign of the total interaction is not the 
same everywhere). Apart from these considerations, it seems worthwhile to examine the electron-phonon 
interaction in the cuprates since there have been some experimental indications fmm Point Contact 
Spectroscopy that the phonons play a crucial role in the superconductivity of YBCOQ. 

We therefore consider a strong electron-phonon interaction through our generalized Eliashberg 
schemeH~U. In this scheme we solve the Eliashberg equationstHi for an electron-phonon interaction hav- 
ing a very strong k-dependence of its matrix elements .g(k, k') that enters in the phonon propagator 
£'(k, k'; — ujn') through eq. (|^). Apart from the angular dependence discussed above, the same argu- 
ments about the short range of the interaction in momentum space also apply to the dependence oapthe 
bare quasiparticles energy (i.e. on the component of the momentum k perpendicular to the FSll3ilil). 



Therefore, we consider an interaction that is cut off for high energy transfers, i.e. ^(k, k') ^ only if 
\£.k — ^k' \ < Cent, where £,cut is the energy of the cutoff. Furthermore, we shall exclude the study of the 
p-wave pairing from our analysis since we shall only consider the case of singlet pairing. 

However, we want to emphasize, that the existence of a d-wave solution of Eliashberg equation does 
not rely on a cutoff at low energyliJ since as we shall see in the following section, even the weak coupling 
treatment in which the ^fc-dependence is completely neglected can exhibit a d-wave solution. Furthermore, 
we show in section IV, that the strong coupling case without cutoff also gives d-wave solutions. 



II. WEAK COUPLING SITUATION 



Although a strong coupling calculation is needed to treat correctly the case of the High- Tc Supercon- 
ductors, the important parameters, as well as the main trend, can be put forward by the standard weak 
coupling theory. We will consider here an isotropic FS, and the electron-phonon interaction as being 
responsible for superconductivity. We want to study the effect of the anisotropy of the interaction, as 
well as of the temperature on the symmetry of the gap in the superconducting state. 

We solved numerically the standard gap equationEI restricted in 2 dimensions 

A{e) = — r de' r"°' d^,,__^E^tanh i Ve + A^(g0 \ ^y^^^r^ _ ^i*{e,e')] (i) 

where ujmax is the maximum phonon frequency. The Coulomb pseudopotential /i* is chosen as con- 
stant for simplicity (in this section). We consider a very anisotropic electron-phonon interaction V{9,6') 
corresponding to a predominantly forward scattering situation. We chose to modelize the behaviour of 
V{9, 9'), where the angles 9 and 9' refer to points k and k' on the FS (cf fig 1), with a gaussian : 

"<".'") ^5J^cxp(-tp!). (2) 

the parameter 59 determining the angular range of the interaction. 

For T = and a small gap A{9) < ijJmax (i-e. a small electron-phonon coupling constant A), eq. (|^) 
can be written 



A(6») = — / d6l'A(e'')ln 
27r 



A{9') 



[V{9,9')~^,*{9,9')]. (3) 



We find that there always exists an s-wave solution to eqs. or but nevertheless a d-wave 
solution can also be found provided the angular range of the interaction 59 is small, i.e. the electron- 
phonon scattering is predominantly forward. We show the shape of the gap lS.d{9) for the d-wave solution 
for 2 different values of 59 in fig. |l]a. We see that the shape of Ad{9) becomes square-like when 59 0, 
whereas it is cos-like for larger 59. The temperature also has the same effect of softening the shape of 
the gap (see fig. |l|b), diminishing thereby the anisotropy. Fig. ||a shows how A^; = maxA(;(6') depends 
on the angular range 59. We see that A^ decreases rapidly with 59, falling to zero when 59 « tt/2. The 
values of Ag for 2 different values of n* are also shown for comparison. We plotted the corresponding 
binding energy Ebmd = ^^2^^ ^ d9A'^{9) in fig. ||b (we set N{eF) — 1 for simphcity). These 2 figures 
show clearly that a large 59 will favour s-wave and a small one the d-wave. The crossover 59c between 
these two symmetries is given by E'f^^^{59c) — e\^'1^^{59 c) . When 59 > 59c, the system will prefer s-wave 
even if the s-gap is smaller then the maximum value of the d-wave gap. Then, in the case presented here, 
the crossover, for fi* — 0.05, will take place at 59c ~ 7r/14, but the maximum d-wave gap will remain 
larger than the s-wave gap until 59 « tt/8. As it can be seen from fig. the temperature dependence 
of the gap is quite BCS-like. We observed that Tc and the shape of A(T) do not depend significantly 
neither on the angular range 59, nor on the symmetry of the gap (this was done by varying T, keeping 
59 constant). Therefore a temperature- induced crossover from d to s-wave, with constant (temperature- 
independent) parameters is unlikely. However, an increase of the temperature might broaden the range 
of the interaction 59 since the predominantly forward scattering may be somehow washed out by thermal 
effects. Then, since a relatively small change in 59 can have a very big effect on the magnitude of the 
gap (as it is seen from fig. ||), the temperature-effects will be seen through the dependence of 59 on T. 
Therefore one has the possibility of a crossover from d to s-wave induced by the temperature. When 59 is 
very small it is clear that other symmetries higher than d can be possible, but here we have not studied 
such solutions. 



We also tried another possibility for the attractive coupling V{9,9'). We considered that the strength 
of the on-site coupling could vary with the angle 9, in order to favour the apparition of nodes in the gap. 
We then chose to replace A in eq. (|) by X(p{9)ip{9'), with (p{9) = cos^" 20 or ip{9) = ^(1 + cos2"26'). 
Indeed the strength of both of these coupling would be smaller at the place on the FS where a node of 
the d-wave solution should be found. The solutions for different sets of parameters showed that although 
the shape of the gap A(0) was altered significantly, the main qualitative results about the existence of a 
d-wave solution remained unchanged, i.e. the important parameter allowing for a d-wave solution is 69. 
Indeed if 69 is too large (> tt/IO), not only the d-wave solution disappears, but the superconductivity as 
well, i.e. A{9) = 0. 

From this weak coupling calculation, we clearly see that the key physical property that allows the 
d-wave symmetry is the tendency to forward scattering. Indeed, since there are only few electrons side- 
scattered, the gap at a point of the FS will only have very little influence on the gap in the neighbouring 
region. Then it is possible to have a very anisotropic solution. 



III. STRONG COUPLING CASE 



The value of 2A/kBTc in the cuprate clearly shows that-4hese materials are in the strong coupling 
regime. We have therefore solved the Eliashberg equational3 to find the superconducting state. Since 
the FS of the cuprates is mainly 2 dimensional, we consider the form of the Eliashberg equations that 
is already integrated along the remaining momentum direction, i.e. these equations specialized to 2 
dimensional momentum space. Instead of using {kx,ky) as momentum variables, we will use the more 
convenient 'jpplaf-|like" k = (0,^fc) where 9 is the angle and = — e^?. As we already did in previous 
publicationsB"llil't2l, we do not integrate out the dependence on the momentum ^fc. The Eliashberg 
equations can then be written in the Matsubara frequencies formalism {uJn — 7rT(2n + 1)): 

n' 

ixiK,u„) - 1) ^ -ty: jj d^/.- ^^"^^^^^"' " "^^^ ""^"^'""'^ (5) 

n' 

[Z[K,UJn) ~l)UJn=T d k' (6) 

n = {Z{K',iJ,,,M^ + {X{K',UJn')^k')^ + $(«',C^n')' (7) 

where D{k^ k' ,lu„ ~ ^n') is the phononic propagator given by (we assume here an Einstein spectrum 
for simplicity) 

£>(/«, k', u;„ - w„/) =5(k,/«')2- '^^''^ — (8) 

and K,') is the Coulomb repulsive pseudopotential that, we will consider as mainly energy- 
independent for simplicity (we think it should depend on ^fc, ^^'cJ; although this effect is not negligible, 
we will neglect it here for simplicity). Therefore, we have chosen /i(K, k') = fi(9^9'). In practice we use 
weighted Matsubara frequencies when solving eq. (^j^) in order to speed up the computationsll3. 

As we already pointed outM, it is possible to explain some typical high properties from a model 
where the electron-phonon interaction is severely cut off along the S^k direction. Indeed, the dielectric 
constant of the ionic background Cion is anomalously big at small energies^, i.e. eioniS, 0) « 40 — 60 
whereas it is normal at higher energies, that is of the order of 2. This will have the effect of reducmg the 
electronic screening at the low energies, thereby enhancing the electron-phonon interaction therdiJ. We 
modelize this behaviour through a cutoff for low energies in the electron-phonon coupling constant: 

g{K, k') = .g(0, 9') e(|a - 6' I - Ut{e, 9' j) (9) 

where Q{x) is the Heaviside step function and ^cMt(^j ^') determines the cutoff and is typically of order 
ujph/2 (except where the elastic scattering rate by defects and impurities r^^^ is big, forcing the cutoff to 
be larger (see below)). 

For the angular dependence of the interaction, guided by the above weak coupling calculation, we 
choose the interaction g{9,9') to be rather strongly peaked around the angle 9, which means that 69 
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in eq. (g) is small. Therefore X{d,6') = 2N{eF)g'^{{d,0')/ujph corresponds to the situation of forward 
scattering that favours anisotropic solutions of the gap. 

In our formulation of the Eliashberg equations, Z{K,uJn) is the standard mass renormalization and 
X{K,ujn) is the renormalization of the bands due to the interaction effects. Because of the small cutoff, 
X depends strongly on the temperature and the gap. The gap is given by 

Z{K.,UJn) 

The strong ^fc-dependence of the interaction has very interesting effects on the gajMiarameter $ and on 
the renormalization functions X and Z that we studied in some length elsewhereBiill. In the context of 
the study of the symmetries of the gap, these effects are not very important because they are somewhat 
hidden by the effects of the anisotropy of the coupling. Nevertheless, we wanted to keep the same 
formalism in order to continue to e}a3iei|re the richness of the solutions it can bring, and in order to be 
consistent with our earlier approachciEJ. 



IV. RESULTS 

We solved the set of equations (|^0) numerically by iterations. For the discretization of the integration 
over ^fe/ we used a mesh with 101 points, and for the integration over 9' , we split our simplified FS (see 
below and fig. ^) in 4 or 8 pieces in an angular range from to tt (one half of the total FS). The suhj. 
over the Matsubara frequencies uJn was carried on a set of frequencies equivalent to the 68 first onesllj 
which means that the range along the imaginary axis waa-at least 8 times the phonon frequency . 

We take the same shape of the FS as in earlier worksM'Ej. This is the simplest case one could imagine 
beyond the very rough approximation of an isotropic FS. It looks like a square with rounded corners (see 
fig. and roajy bear some resemblance with the FS observed by ARPES measurements in YBCOEa 
or in BiSCCOEj. We consider a FS consisting mainly of 2 different kinds of pieces: the planes and the 
corners. The corners are near a Van Hove singularity (thereby having a large density of states) and vp 
(perpendicular to the z-axis) is then much smaller than on the planes. Consequently the elastic scattering 
rate is larger on the corners, smearing out the electronic momentum (^fe) and preventing the cutoff 
icut from being very small. This is why we will choose a much bigger cutoff for the corners than for the 
planes. In order to be able to reproduce all usual symmetries of the gap (i.e. s and d-wave), we need to 
split the FS in at least 4 pieces in a range of tt (half the FS) . Fig. ^ shows the splitting of the FS in 8 
pieces, the 4 pieces case being simply obtained by grouping the 3 pieces on the "plane" part of the FS in 
a single piece. 

For simplicity we set LUph — 1 and N{k) = 1 V /« (the K-dependence of N{k) is implicitly incorporated 
in g{K, k')). All the energies are given in units of coph. Since we only consider a small number of pieces, 
specified by their angular position 9i,i = 1, ... ,4 or 8, we will use a matrix notation, denoting g'^{Oi, 9j) 
by gfj. All these matrices gfj,^,cutij and /iy are obviously symmetric. We start with the description of 
the parameters of the 8 pieces case, the case with 4 pieces follows the same general philosophy and its 
parameters will be given subsequently. As was mentioned above, the cutoff ^cut on the corners cannot 
be as small as the one on the planes. We have chosen ^cutn = Cc«t55 = 2, and £,cutii — 0.4 for all other 
diagonal matrix elements, i.e. for i 1 and 5; the non-diagonal elements were set to —2 \/i ^ j. 
The angular dependence of the electron-phonon coupling = 2gfj (for our choice of uiph and N{k)) is 
basically as described by eq. (||). The situation is just a little bit complicated by the anisotropy of the 
FS and our choice for the splitting in which the pieces have not the same weight. The diagonal elements 
were chosen as follows: An = A55 = 2 (corners), A22 = A44 — Aee = Ass = 1-2 (part of the planes 
neighbouring the corners), and A33 = A77 = 0.6 (central part of the planes). The non-diagonal elements 
are rapidly decreasing with the distance in momentum space: A12 = Ais = A45 = A56 = 0.7 for the 
coupling between neighbouring corners and planes, A23 = A34 — Ag? — Xrs — 0.4 between neighbouring 
parts of the planes, and A^ = 0.2 for the rest. This corresponds approximately to a gaussian shape in 
the angular interaction with 59 ~ tt/8. The Coulomb repulsion is, as we already said, mainly constant, 
but with only some variations taking in account the different weights of ithe different parts of the FS and 
the repulsive interaction between corners as seen in neutron experiment^3and described by the theory of 
Monthoux et aU. Therefore, we have a more repulsive interaction (that will favour the d-wave solution) 
between electrons belonging to pieces separated by an angle A9 = \9 — 9'\ of about 7r/2. In the light of this, 
we set the matrix elements of fj, as follows: the corner-corner interaction fj,c-c{AS = '^/'^) = Mis = 0.4, 
the plane-plane one fj,p-p{A9 = 7r/2) = fi26 — ^J■37 — fJ-is — 0.3, the corner-plane one (for pieces separated 
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by 37r/8) fic-p{'^d ~ Svr/S) = /X14 = /xig = /U25 = 0.3, and the rest Hij = 0.1 except for the 2 diagonal 
elements standing for the on-site repulsion of the corners (Xc-c{^d = 0) = /in = /i55 = 0.25 (to take in 
account the weight of these pieces) . 

The results of this calculation with 8 FS pieces are shown in fig. ^ for T ^0 case (for technical reasons 
the temperature was set very low, T = 0.02wpft). The gap function <&(/«, ujn), as well as the renormalization 
functions w„) and Z{k, cun) are plotted as functions of for each of the 8 FS pieces (denoted by 
the angle 8j = ^(j — 1), j = J- . . . , 8) and for n — 0, since we believe that the low to behaviour of these 
functions is roughly constantO, so that it is satisfactory to have their values at w = (which in turn is 
approximately the value for the w„=o — Matsubara frequency). Both s- and d-wave solutions were 
found for this set of parameters (dashed and continuous lines respectively), the symmetry of the solution 
depending only on the choice for the symmetry of the initial condition. The true thermodynamic gap of 
the s-wave solution — $/Z (eq. ([l^)) is quite anisotropic, although less than its function <I>, since 
the behaviour of the function Z tends to diminish this anisotropyQ. There are 3 different values for the 
gap As at the FS: A^i = A^s = 0.48, As2 = A34 = A^e = A^s = 0.35 and A33 = A^y = 0.32 (in units 
of ujph), where A^j = As{9j^^k = 0,a;„=o)- On the other hand, the maximum value of the gap for the 
d-wave solution A™"^ is much larger, i.e. A™°^ = ma.x Ad{0j) = 0.88 (on the corners). We could then 
expect from our weak coupling considerations that the system will prefer the d-wave symmetry since such 
a large difference between the magnitudes of the s- and d-wave gaps indicates that we are not so close to 
the crossover described in fig. ||. 

In order to show the effect of the choice of the initial condition, we studied a case with only 4 FS 
pieces (grouping the 3 pieces of the plane in 1) for 3 different initial conditions. We considered an s-wave 
starting gap $o,s('*i ^n) — ^ V K,n and 2 d-wave ones corresponding to a dx2-y2 , i.e. $0.^ 2_ 2 1**! ^n) = 
cos20 V ^k,n- and to a d™, i.e. $o^d^^(K,ti;„) — sin2^ V ^fc,n. The solutions for the gap parameter 
$(/«,a;„) are shown in fig. o. The coupling parameters are set with the same philosophy as for the 8 FS 
pieces case and are as follows: 

/ 0.25 0.25 0.4 0.25 \ 
0.25 0.25 0.25 0.4 | 
0.4 0.25 0.25 0.25 ' 
\ 0.25 0.4 0.25 0.25 / 

where the pieces numbered 1 and 3 stand for the corners and the two others for the planes. Here we 
chose a somewhat less anisotropic electron-phonon interaction than for the 8 FS pieces case. The range 
of the interaction can be estimated as 69 ~ n/S (but the angular dependence is not exactly gaussian). 

The solution shown in fig. ^ possesses an s-wave symmetry, as expected. The anisotropy of the gap is 
similar to the one of the calculation with 8 pieces. The true gap A^ as defined in eq. ( p^ ) is (for ^fe = 
and n = 0) Agi = A^a = 0.44 on the corners and As2 = As4 = 0.32 on the planes. The d^2_y2 solution 
in fig. ^ also has the same symmetry as the initial gap. In fact, this solution has not exactly the d,j.2_y2 
symmetry since the nodes are not found at 9 = tt/A but they are probably slightly displaced somewhere 
in the seams between pieces 1 and 2, and 3 and 4 (corner and plane). This is due to the small number of 
pieces considered. The maximum value of the gap is smaller than for the corresponding s-wave solution, 
i.e. AJ^"^^ = 0.36, as it can be expected since the angular range of the interaction 69 is larger. The last 
result presented in fig. |^c is very similar to the previous one. The only visible difference with the gap 
parameter of fig. ^ is that the sign of its component on the plane (which is small) is reversed. This 
seems to be the only effect of the dj^y-like initial condition. 

We have also observed asymmetrical solutions of the Eliashberg equations (neither d nor s-like), i.e. 
of lower symmetry than the interaction and the FS that indicate a spontaneous symmetry breaking (a 
kind of Jahn- Teller effect). We will discuss the physical meaning of these solutions in a forthcoming 
communication . 

Finally we also considered the 4 pieces without cutoff parameters in order to see if the question of 
s- vs d-wave symmetry was affected by this. The new solutions are qualitatively the same as when the 
cutoff was introduced. Clear d-symmetry are found, but there is no dispersion of (f> with ^k- The plots 
corresponding to fig. |^ consist just of straight, horizontal lines, and we do not show them here. Thus we 
conclude that the cutoff parameters are not crucial for the appearance of d-wave vs s-waves. Still, the 
cutoffs lead to other characteristics of the gap etc. 
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V. DISCUSSION 



As we already noted in a previous workEj, the choice of the symmetry of the initial gap used to start 
the iterative procedure is crucial for the symmetry of the solution. Indeed, there always exists an s-wave 
solution of the Eliashberg equations (^||), but a c?-wave solution can also be found (depending this time 
on the values of the other parameters) provided that the iterations were started with a rf-wave like gap. 
This duality of solutions can be easily understood. In the limiting case where SO ^ (and /i — > 0), i.e. 
where each FS piece is entirely decoupled from the others, the phase of the gap parameter $ can be chosen 
freely on each piece. Then there is no constraint on the symmetry of the gap, or, in other words, the 
s-wave and d-wave solutions (or sometimes even higher symmetry) are completely degenerate. In the case 
considered here, the interaction between electrons of different pieces, although not strictly zero, is quite 
small. We are then in a nearly degenerate state where both s and d wave solutions can coexist, but where 
the relative phase of the gap parameter between pieces is determined. In the light of this discussion, it 
seems rather astonishing that we were not able to get a clear c?2,j^-wave solution (fig. ^c). Indeed, the 
matrix A^- is symmetrical with respect to the exchange of plane and corner, and therefore the symmetry 
of the result should depend on the symmetry of the initial condition. But the asymmetry between planes 
and corners is in the cutoff matrix , reducing sufficiently the average attractive coupling on the planes 
to impose the dj,2_j,2-wave solution over the d^y one. (A similar effect can be obtained without cutoff by 
setting smaller values for the diagonal element Xa of the planes than for the one of the corners). 

The angular dependence of the CoiJ£)mb pseudopotential fj,{6,6') was set in order to mimic the be- 
haviour observed by neutron scatteringc2l and the theory Monthoux et al.El. This favours a d-wave solution. 
However, this enhancement of the repulsive interaction for \9 — 6'\ « 7r/2 is not essential for obtaining 
a d-wave. Indeed other calculations in the 8 FS pieces case with a constant Coulomb pseudopotential 
Mij = Mo ^hj exhibit a d-wave solution quite similar to the one presented in fig. ^ From another 
calculation (leading to an s-wave solution), it seems that the condition for having a d-wave solution is 
that Ay < 2/iy when \9i — 0j \ ~ tt/2. 

The effect of the repulsive potential /i on the symmetry of the solution can be seen clearly by comparing 
fig. ||a and ^ (or the two solutions on fig. ||). It can be noted that under the action of the repulsion, 
the gap parameter $ tends to have nodes, i.e. to have different signs at different places of K-space. In 
the d-wave case (fig. ||b), the nodes lie on the FS, located at 9i « f (2n -I- 1), and the sign of the gap is 
mainly constant along (for the FS pieces very near the nodes, there are other nodes removed from the 
FS due to fluctuations and the small value of the gap on these pieces). On the contrary, in the s-wave 
case (fig. ^), the sign of the gap does not depend on 9 and the nodes are removed from the FS, the 
change of sign of the gap happening along the S.k direction. 

In what concerns the nodes for the d-wave solution, we want to emphasize that they are not exactly 
located at 9 ~ tt/A (see fig. ^ and ||). This is due to the relatively small number of FS pieces considered in 
our calculations. But iUis nevertheless interesting to notice that this reproduces quite well the situation 
observed by Ding et aO in BiSCCO. Indeed they observe the nodes at 9 = j±6ip and this is compatible 
with a situation where the experimental result is the average of two degenerate d-wave states, one with the 
nodes shifted by Sip and the other by —Sip. Such a degeneracy may be due to different orientation of the 
different grains or to orthorhombic distortions. However, the d-wave symmetry for these materials is still 
subject to controversy and this explanation does of course not exclude the possibility of an anisotropic 
s-wave or of some effect of the lattice modulation in BiSCCO. 

A salient feature of the present model is, that when the d-wave superconductivity is destroyed, for 
example by defects or impurities that scatter the conduction electrons by a large angle, then s-wave 
superconductivity takes over, without a large change in the superconducting gap, or in Tc- Experimentally, 
we do not know yet whether, for example, YBCO with defects is an s-wave superconductor; but some 
experimental^ favour such a picture. 

In this picture, the Van-Hove (or extended Van-Hove) singularity (VHS) near the FS is a secondary 
factor. The d-wave pairing exists even for an isotropic FS, without any VHS, as long as the|-scattering 
angle 69 is sufficiently small (section ||). The small S9 also accounts for the high value of Tcllil. In this 
respect, our model differs-i'rom the theory of Pines and collaborators, and the more recent extension of 
this theory by AbrikosovcHl. When a VHS is present near the FS, it favours the d-wave pairing, causing 
it to prevail at a somewhat larger value of 69, and it further increases the value of Tc. We believe that 
the increase in Tc of HgBaCaCuO from 145 K at ambient pressure, to 164 K under pressure observed by 
Chu et aS, is due to the VHS approaching the FS. However, the high value of Tc = 145 K at ambient 
pressure is due predominantly to the small value of 69. 
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VI. EXPERIMENTAL SUPPORT OF THE PROPOSED MODEL 



The main assumption of this work is, that the superconducting pairing is caused by phonons, and 
that these phonons have a very small q value, q/2kp sa tt/IO. Several experiments lend support to this 
assumption. 

1. A structure is found at anj-energy of about 40 meV above the gap in YBCoi and about 45 meV 
above the gap in BiSCCOcJ. Vedeneev et al interpret it as a McMillan-Rowell structure due to 
phonons, and extract from it a value of A consistent with the value of Tc- 

2. ThCjiialue of A extracted from the resistivity of YBCO (and other cuprates) is very small Martin 
et a£5 find A « 0.1. The value of A found from band-structure calculations is about 1.7l3 and the 
nmaimum value needed to account for Tc of YBCO (assuming a very small value of fj,*) is about 
3cZl. The transport A is a factor of « 69"^ smaller than the McMillan A (this is essentially the 
1 — cos6' factor entering the "classical" expression for the resistivity). This indicates that if the 
superconductivity is indeed due to phonons, S9 must be rather small. 

3. A charse. density wave (CDW) with wavevector q = 0.24 A"-'^ is observed in BiSCCO overdoped with 
oxygenEa. Since kp « 0.7 A-\ q/kp « 0.35. This value of q is close to the value of the Thomas- 
Fermi screening parameter, ktf = i^tf/ y/^ion- Thus, 69 sa tt/IO. This is close to the value of the 
crossover from d-wave to s-wave pairing described in section O. Obviously, the estimates made in 
the present work are rather crude and not quantitative. Nevertheless, the small value of q provides 
strong support for the present model. 

We wish to point out, that while the observation of a CDW is a proof that the electron-phonon 
interaction is very strong, the converse is not necessarily true; a very strong electron-phonon interaction 
causes an appreciable softening of the phonons, but not necessarily to w = 0, required for a static 
distortion. Therefore, the absence of the CDW in optimally-doped BiSCCO does not indicate that 
the electron-phonon interaction there, with a, q ~ 0.24 A^^ phonon, is not very strong. One possible 
mechanism to account for the difference in behaviour between optimally-doped and oxygen-overdoped 
BiSCCO, is as follows: In optimally-doped BiSCCO, the strong electron-phonon coupling softens the 
LA phonon (for values of q around 0.24 A~^), but not to zero, thus there is no static distortion. The 
inherently soft TO phonons (with frequencies around 10-15 meV) possess a different symmetry, and 
are thus not strongly affected. In oxygen-overdoped BiSCCO, the interstitial oxygens in the (Cu02)n 
planes destroy the symmetry, and consequently the LA and TO phonon modes are coupled, and one such 
coupled mode softens to zero, producing the static distortion (CDW). Thus, the high Tc, and pairing 
with d-wave symmetry, caused by the electron-phonon coupling of the LA phonon with small q, do not 
necessitate a static distortion; but the static distortion of the overdoped material points to the very strong 
electron-phonon coupling at this q value. 

VII. CONCLUSION 

We have shown in this paper that a phonon-mediated mechanism can lead to a solution of the Eliashberg 
equations with d-wave symmetry, contradicting thereby the conventional belief that it can only lead to 
s-wave solutions. This belief is based on the fact that the electron-ion interaction is nearly isotropic in 
normal metals. However the situation is quite different in the cuprates since the large polarizability of the 
oxygen ionic background prevents an effective electronic screening. Therefore, the range of the interaction 
in momentum space is much smaller than in normal metals. This is sufficient to allow for d-wave solutions 
in the weak coupling case (the Coulomb pseudopotential fj,* being constant, it gives no contribution to 
the integral in eq. (|l|) for a d-wave gap). However, the strong coupling case requires an appropriate 
repulsion, i.e. a total repulsive interaction, for |0 — 0'| « tt/2 in order to allow the d-wave symmetry. This 
extreme sensitivity to the parameters shows that our solutions of the Eliashberg equations are nearly 
degenerate and that a small change in the coupling parameters can cause a crossover between the two 
symmetries. 
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FIG. 1. Shape of the gap Ad{0) with the couphng parameters A = 0.5 and /i* = 0.02. (a) For 2 different 
angular ranges 59 = 7r/30 (plain line), 59 = tt/IO (dashed line) and a cosine function is also plotted for comparison 
(pointed line); and (b) for different temperatures, i.e. T = 0.01, 0.09, 0.11, 0.13, 0.15 and 0.29 (with 59 = tt/SO). 
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FIG. 2. Dependence of the solution on the anguhir range 56. (a) Maximum of the rf-wave gap A^i vs 56 and 
values of the s-wave gap for ji* = 0.05 (dashed line) and /x* = 0.15 (dash-dot line), (b) Binding energy (with 
N(eF) = 1) for the same cases as before. 
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FIG. 3. Temperature-dependence of the maximum of the gap for the rf-wave case with 69 = n/30 (plain hne 
and circles) and for the s-wave with fj,* = 0.02 (dashed line and crosses). 
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FIG. 4. Schematic view of the FS of the cuprates. It looks more or less like the one observed in the cuprates. 
We also show the labelling of each FS piece in the case of 8 pieces. 
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FIG. 5. Solutions of the Eliashbcrg equations $, X and Z as functions of the bare quasiparticles energy for 
each of the 8 FS pieces. The s-wave symmetry (dashed hne) and d-wave symmetry (plain line) are shown. The 
parameters describing the electron-phonon coupling and the Coulomb pseudopotential are discussed in the text. 
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FIG. 6. Solutions of the Eliashbcrg equations $ as function of the bare quasiparticles energy for each of 
the 4 FS pieces and for LUn=o- The 2 pieces numbered 1 and 3 correspond to the corners and the 2 others to the 
planes. Different initial conditions for the gap are considered, (a) With an s-wave initial gap, (b) with a d^2_y2 
one, and (c) with a dxy one. The temperature was very low, i.e. T = 0.02 and the coupling parameters are given 
in the text. (All energies are in units of uiph)- 
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